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Abstract 

Gluonic correlation functions and Bethe Salpeter amplitudes are calculated 

in an instanton-based model of the QCD vacuum. We consider both the pure 

gauge case and the situation for real QCD with two light quark flavors. We 

show that instantons lead to a strong attractive force in the J^*^ = 0^^ 

channel, which results in the scalar glueball being much smaller than other 

glueballs. In the ^ channel the corresponding force is repulsive, and in 

the 2"*"^ case it is absent. Due to the strong classical field of the instantons, 

these forces are even stronger as compared to the ones for mesons made of 

light quarks. The resulting correlators, masses, coupling constants and wave 

functions are compared to results of lattice simulations. 
PACS numbers: 12.38.Lg, 12.38.Gc, 14.40.Cs. 
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One of the most obvious question in QCD is why the observed hadrons are made of 
quarks and not of gluons. It appears that glueballs are much heavier than typical quark- 
model hadrons, and therefore they have large widths and/or complicated decay patterns, 
making them difficult to find. But why are glueballs so heavy? What are their masses, radii 
and other parameters in a purely gluonic world, and how do they change if one includes 
light quarks? 

To get a reference point, consider a conventional model of glueball states such as the bag 
model. The lowest fermion and (electric) gluon modes in a spherical cavity have energies 
2.04/i? and 2.7/i?. Thus, glueballs are expected to be heavier than quark states, but not 
much. Ignoring spin-dependent forces, one expects mo++ — m2++ — 1 GeV and mo-+ — 1.3 
GeV. Including these forces (and other refinements) [|l|] the model predicts that the low-lying 
glueballs have masses m ~ (1.0-1.8) GeV and very similar radii r ~ (0.7-0.9) fm. 

A number of "glueball candidates" have been experimentally observed, but none was un- 
ambiguously identified (see however [0). Lattice simulations provide important qualitative 
insights, and (although large-scale numerical efforts are still necessary) a few statements 
appear to be firmly established 0]: (i) The lightest glueball is the scalar, with a mass in 
the 1.6-1.8 GeV range, (ii) The tensor glueball is significantly heavier m2++/mo++ — 1.4, 
with the pseudoscalar one heavier still mo-+/mo++ = 1.5-1.8 |Q. (iii) The scalar has a much 
smaller size than other glueballs. This is seen both from the magnitude of finite size effects 
and directly from measurements of the wave functions The size of the scalar glue- 

ball (defined through the exponential decay of the wave function) is ro++ — 0.2 fm, while 
r2++ — 0.8 fm 0. For comparison, a similar measurement for the vr and p mesons gives 0.32 
fm and 0.45 fm |^, indicating that spin-dependent forces between gluons are stronger than 
between quarks. 

Important tools that provide information about gluonic interactions are the correlation 
functions of gluonic operators with the relevant quantum numbers, such as the field strength 
squared (0++), the topological charge density (0"+), and the energy density (2++): 
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Os — (gG'^^u)'^, Op — -^fiupag'^G^^u^pa, 
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Ot = -^{gG'^^uf - g^Gi^Gi^ . (1) 



In the following, we consider correlation functions nr(x) = (0|Or(x)Or(0)|0) for euclidean 
separation x. An important low energy theorem was proven in the integral of the 
scalar correlator is determined by the gluon condensate, / d^x 115(0;) = ^^^""^ {(gCy), where b 
denotes the first coefficient of the beta function and the integral is regularized by subtracting 
the perturbative contribution. This theorem indicates the presence of rather large non- 
perturbative corrections in the scalar channel. On the other hand, the operator product 
expansion (OPE) predicts that the leading-order power correction 0{{G'j^J) /x^) vanishes 0, 
while radiative corrections of the form as log(x^) {G"^^^) /x^^ or higher order power corrections 
like {gr^'^Gl^GlpG''p^^)/x^ are very small. 

In practice, there are two approaches to QCD sum rules for scalar glueballs. In the 
original work the low energy theorem was enforced by introducing a subtraction constant. 
In this case, the subtraction constant completely dominates over ordinary power corrections, 
and one finds a large scalar glueball coupling. However, the vacuum picture advocated in this 
paper was a rather homogeneous one (with large size instantons melted with other vacuum 
fiuctuations), so that the source of the large subtraction constant was not clear. In more 



recent works on the subject [|,1C], the low energy theorem was not enforced and instead a 
number of higher order corrections in the OPE were evaluated. Although the resulting mass 
estimate is similar to what was obtained earlier, the resulting glueball correlation functions 
are very different. In particular, the low energy theorem is underestimated by about an 
order of magnitude and the scalar glueball coupling constant is substantially smaller. 

In this paper we study gluonic correlation functions in an instanton based model of the 
QCD vacuum [|TT|]. The main assumption underlying the model is that the gauge fields 
in the QCD vacuum are dominated by the strong fields of small size instantons. Support 
for this assumption is provided by the calculation of hadronic correlators in the model ||T^ 



and the analysis of "cooled" lattice configurations ||13|, where non-classical fiuctuations are 



removed from the vacuum. Our main point is that small size instantons can lead to a strong 
enhancement in the scalar correlation function and give a consistent description of the low 
energy theorem. We believe that the smallness of the scalar glueball provides a strong 
argument in favor of the picture presented here. 

Before we discuss quantitative predictions, let us consider qualitative features of point-to- 
point correlators . Quark correlation functions roughly fall into three classes, depend- 



ing on how asymptotic freedom is broken at intermediate distances. The non-perturbative 
corrections may either be (i) large and attractive (the corrections have the same sign as the 
free correlation function), as for the pseudoscalar vr, a mesons; (ii) large and repulsive, as 
for the heavy scalars r( and 5; (iii) or they can be small even at rather large distances x ~ 1 
fm, as is the case for the vector mesons p, 01,0;, 0. This classification is easily understood 



since the instanton-induced interaction between quarks found by 't Hooft has precisely 
the required spin-isospin properties. 

Instanton effects in the gluonic correlation functions can be studied by calculating the 
correlator in the classical field of a single instanton. Adding the short distance contribution 
from the free gluon propagator, one finds 
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n5,p(x) = (±)^^ + np"ni„,t(x), (2) 
nT a: = (3) 



where g is the running coupling constant and Tlinstix) is the instanton contribution 

, , 122887r r s fi A 9 
n.„.,(x) = y,^y,^\^A y + - 9V - 160y^ - 120 

with y = x/ p. The approximation used is that we ignore any interference between quantum 
and classical fields. The reason is that these contributions correspond to power corrections 
0((G'^j,)/x^), which, as mentioned above, are very small in the channels considered here. 

To first order in the instanton density, we find the three scenarios discussed above: 
attraction in the scalar channel, repulsion in the pseudoscalar and no effect in the tensor 



channel. The last case is a consequence of the fact that the stress tensor in the self-dual 
field of an instanton is zero. If one compares the result (4) with a similar calculation for 
the pion, one finds that the instanton contribution in the glueball channel is enhanced by 
a factor Sq = {Sn^ / g{pY)^ with respect to the result for the pion. This means that despite 
the fact that the scalar glueball is so much heavier than the pion, the correlation function 
at distances x ~ p is even larger. 

In order to make these statements more quantitative and study their dependence on 
the presence of light quarks, we have calculated the correlators for three different instanton 
ensembles. The simplest is the Random Instanton Liquid Model (RILM), which assumes 
that instantons and anti-instantons are distributed randomly in position and color space. 
Already this simple model is very successful in the description of a large number of hadronic 



correlation functions In this model the QCD vacuum is dominated by small-size (p ~ 0.3 
fm) instanton or anti-instanton fluctuations with a total density n ~ 1 fm~^. Including the 
correlation between instantons introduced by the gluonic interaction between them, gives a 
more complicated ensemble which we call the Quenched Instanton Liquid Model (QILM). 
Also taking into account the fermionic determinant one arrives at the unquenched Interacting 
Instanton Liquid Model (IILM). As shown in |T^, this ensemble reproduces an important 



feature of QCD : the screening of the topological charge and a correct description of the rj' 
channel. 

In our simulations we calculate the correlators as in (H), but the classical part is now 
evaluated for multi-instanton configurations. The scalar correlator has an x-independent 
contribution from the gluon condensate, which must be subtracted. The running coupling 
constant g is calculated from the perturbative beta function at short distances, but frozen 
at a value of a^/vr = 0.3. 

The resulting correlation functions are shown in fig.l. The correlators are normalized 
to the free ones, so that all curves approach one at short distances. Deviations from one 
at intermediate distances are very different in different channels. Up to x ~ 0.25 fm, these 
deviations are consistent with the single- instanton correction (Q), but at larger x multi- 



instanton contributions become important. Note that for the pseudoscalar correlator in 
the interacting ensemble the correction even changes sign. This is a result of correlations 
between instantons and anti-instantons that lead to the screening of the topological charge. 
In the ensemble with light quarks a new state, the 77' appears. Apart from this, one observes 
that the results for the three ensembles considered are rather close, suggesting that single- 
instanton effects (rather than correlations among them) are dominant. 

One consistency check is provided by the low energy theorem. We find that the integral 
of the scalar correlation function (integrated up to X = 0.7 fm) is 97 ± 6 GeV^ for the 
random and 66 ± 7 GeV^ for the interacting ensembles, to be compared with the low energy 
theorem value 62GeV^. In fig.l we also compare our results with predictions from QCD 
sum rules. The dotted lines correspond to the glueball parameters obtained in (the 
results from |^ are very similar). We clearly observe that QCD sum rules do not predict a 
substantial enhancement of the scalar correlator, the value of the sum rule is only 13 GeV^. 
Unfortunately, most lattice simulations use very non-local operators (in order to increase 
the ground state signal) and their correlators cannot be compared directly to our results. 
We strongly encourage direct measurements of point-to-point correlation functions and the 
corresponding coupling constants, similar to the results for quark correlators reported in 



We have fitted the glueball masses and coupling constants using a parametrization of the 
spectral function that consists of a zero-width pole and a continuum starting at a threshold 
Sq. For r = S", P the correlator reads 



where D{m,x) = 'm/{AiT'^x)Ki{mx) is the euclidean scalar propagator and Ar = (0|Or|0''^'") 
is the coupling of the resonance to the gluonic current. Statistical fluctuations at large 
distances limit the accuracy of our mass determination, but the coupling constants are 



gluonium in the random model we find a mass mo++ = 1.4 ± 0.2 GeV with the coupling 



0- 




(5) 



determined rather well. Fitting the parametrization (^ to the measured correlator for scalar 
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strength Ao++ = 17.2 ± 0.5 GeV^. In the quenched and the full ensemble the correlation 
function is somewhat smaller at intermediate distances, a consequence of the low energy 
theorem discussed above. At distances x > 0.5 fm there are large uncertainties due to the 
subtraction. The mass and coupling constant are mo++ = 1.25 GeV and Ao++ = 15.6 GeV^ 
in the full theory and mo++ = 1.75 GeV and Ao++ = 16.5 GeV^ in the quenched case. These 
values are about twice as big as the value obtained in the only lattice measurement of this 
quantity, Ao++ = 7.8 GeV'^ ]I8|] (which, according to the authors, is compatible with the low 
energy theorem for glueball masses m ~ 1 GeV on the low end of their statistical errors). 

In the pseudoscalar case the classical and one-loop contributions have opposite signs. 
At distances where they tend to cancel each other our approximation (which neglects the 
interference between the two) becomes questionable. However, the rapid downturn directly 
translates into the position of the perturbative threshold, for which we find y/so ^ 3.0 GeV 
in the random model and ~ 2.4 GeV in the quenched theory. We see no clear evidence 
for a pseudoscalar glueball state below the continuum threshold. In the unquenched theory 
we observe the rj' signal with m^/ ^ 800 MeV and A,,' ^ 7.0 GeV^. Using the anomaly 
equation, this corresponds to = 200 MeV. 

The tensor channel has no large non-perturbative corrections, because isolated instantons 
and anti-instantons have a vanishing energy momentum tensor. Thus the classical contri- 
bution is entirely due to the interaction between instantons. One may therefore question 
the importance of instantons in this channel. On the other hand, the OPE also does not 
predict any power corrections at 0{G^) and 0{G^), and in a self dual background field all 
power corrections vanish the smallness of the non-perturbative corrections may therefore 
survive. We find a small classical contribution to the tensor correlator (see fig.l), but have 
not made an attempt to determine the corresponding mass value. The data can be used to 
put an upper limit At ^ 0.6 GeV^ on the tensor coupling. 

Since the scalar correlator is so much bigger than the other ones, one may speculate that 
the scalar glueball should also be much more compact. We have checked this statement by 
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calculating the Bethe Salpeter amplitudes (or 'wave functions'), defined as 

^^^^^ = il^ n^^o|^r(x)OF(o)|o), (6) 

where O|(0) = (7^G'°^(— ?//2)G^j,(|//2) etc. are point-split operators and y is orthogonal to 
X. By definition, ipi^]) = 1. At small separation x the wave function of the scalar glueball 
essentially measures the size of the instanton, ipsiu) = 1 ~ (2l//3p)^ + 0{y*). 

Our results for the random ensemble are shown in fig. 2 (those for other ensembles are very 
similar). The scalar wave function is indeed found to be very compact. It is not exponential 
at short distances (presumably due to the lack of short range perturbative interactions), but 
the overall shape can be described by an exponential decay ip{y) = exp{—y/R), with a fitted 
radius R = 0.21 fm. This value is in good agreement with the lattice result R ~ 0.2 fm [0. 
The tensor wave function is much larger in size, R = 0.61 fm, to be compared with the lattice 



result R ~ 0.8 fm. Together with our earlier work |T2|, in which we determined the sizes of 
the pion r = 0.56 fm and rho meson r = 0.70 fm, this shows that the instanton model leads 
to significantly larger spin-splittings for the glueball radii. For the pseudoscalar glueball 
the interaction is repulsive and we do not find a localized wave function in our model. We 
therefore conjecture that lattice measurements should find a dip in the wave function at 
small distances. Let us conclude by noting that the observed hierarchy of sizes, from a very 
small scalar to a large tensor and a presumably large pseudoscalar, is of great significance 
for phenomenological searches, since it may affect the branching ratios into different final 
states (see p| and the discussion in [0). 

In summary, we have shown that instanton-induced forces between gluons lead to strong 
attraction between gluons in the O^"*" channel, strong repulsion in the ^ channel and no 
short- distance effects in the 2++ channel. We have calculated point-to-point correlation 
functions using the "instanton liquid" model (with and without quarks). The fitted masses 
are compatible with lattice results. More importantly, our large scalar coupling constants 
are in agreement with the low energy theorem, and the scalar gluonium size is as small as 
0.2 fm, as it was observed in 0]. We have argued that measurements of gluonic correlators 
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and wave functions provide important insights into the structure of the QCD vacuum, and 
lattice measurements of the correlators predicted in this paper would be very desirable. 

We would like to thank P. van Baal for a useful discussion which triggered this paper. 
This work was supported in part by the US DOE grant No. DE-FG02-88ER40388 and the 
A. V. Humboldt foundation. 
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FIGURES 



FIG. 1. Scalar, pseudoscalar and tensor glueball correlation functions normalized to the 
corresponding free correlators. The results in the random, quenched and full ensembles are denoted 
by stars, open triangles and solid squares, respectively. The solid lines show the parametrization 
described in the text, the dashed line the dilute instanton gas approximation, and the dotted line 
the QCD sum rule calculation [10]. The horizontal line in the second figure was added to guide 
the eye, the vertical scale in the third figure is 10~^. 

FIG. 2. Scalar and tensor glueball Bethe Salpeter wave functions in the random instanton 
ensemble. All wavefunctions are normalized to one at the origin. The solid and dashed lines show 
a parametrization of the data used to extract the mean square radii while the dashed curve shows 
the scalar wave function in the dilute instanton gas. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9410372v2 
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